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Abstract 

^— » . 

,£3 I We investigate lattice simulations of scalar and nonabelian gauge 

p ■ fields in Minkowski space-time. For SU (2) gauge-theory expectation 

<D . values of link variables in 3+1 dimensions are constructed by 

'"^ \ a stochastic process in an additional (5th) "Langevin-time" . A 

sufficiently small Langevin step size and the use of a tilted real- 

^ . time contour leads to converging results in general. All fixed point 

solutions are shown to fulfil the infinite hierarchy of Dyson-Schwinger 
identities, however, they are not unique without further constraints. 
For the nonabelian gauge theory the thermal equilibrium fixed point 
is only approached at intermediate Langevin-times. It becomes 
more stable if the complex time path is deformed towards Euclidean 
space-time. We analyze this behavior further using the real-time 
evolution of a quantum anharmonic oscillator, which is alternatively 
solved by diagonalizing its Hamiltonian. Without further optimization 
stochastic quantization can give accurate descriptions if the real-time 
extend of the lattice is small on the scale of the inverse temperature. 



1 Introduction 



Lattice gauge theory provides important insight into strongly interacting 
theories such as quantum chromodynamics (QCD). Typically, calculations 
are based on a Euclidean formulation, where the time variable is analytically 
continued to imaginary values. By this the quantum theory is mapped 
onto a statistical mechanics problem, which can be simulated by importance 
sampling techniques. Recovering real-time properties from the Euclidean 
formulation is a formidable problem that is still in its infancies. Direct 
simulations in Minkowski space-time would be a breakthrough in our efforts 
to resolve pressing questions, such as early thermalization or the origin of 
seemingly perfect fluidity in a QCD plasma at RHIC £Q. For real times 
standard importance sampling is not possible because of a non-positive 
definite probability measure. Efforts to circumvent this problem include 
mimicking the real-time dynamics by computer-time evolution in Euclidean 
lattice simulations |21E1- A problem in this case is to calibrate the computer 
time independently of the algorithm. Another procedure amounts to separate 
a positive factor from the Boltzmann factor to be used for importance 
sampling, and re-weight the configurations with the remaining (complex) 
part. These so-called "reweighting methods" suffer from the problem of 
large cancellations of contributions, induced by the oscillating signs in the 
weight ( "sign problem" ) and from the difficulty of ensuring a sufficient overlap 
between the simulated and the target ensemble ("overlap problem"). 

Direct simulations in Minkowski space-time, however, may be obtained 
using stochastic quantization techniques, which are not based on a probability 
interpretation jU El- In Ref. [S| this has been recently used to explore 
nonequilibrium dynamics of an interacting scalar quantum field theory. 

In real-time stochastic quantization the quantum ensemble is constructed 
by a stochastic process in an additional "Langevin-time" using the 
reformulation for the Minkowskian path integral |HJ E]: The quantum 
fields are defined on a <i-dimensional physical space-time lattice, and the 
updating employs a Langevin equation with a complex driving force in an 
additional, unphysical "time" direction. This procedure does not involve 
reweighting, nor redefinition of the Minkowski dynamics in terms of an 
associated Euclidean one. Though more or less formal proofs of equivalence 
of the stochastic approach and the path integral formulation have been 
given for Minkowski space-time, not much is known about the general 
convergence properties and its reliability beyond free-field theory or simple 
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toy models [HI E]. More advanced applications concern simulations in 
Euclidean space-time with non-real actions ^01 [n]. Besides successful 
applications, major reported problems in this case concern unstable Langevin 
dynamics and incidences of convergence to " unphysical" results ^21 EE E] • 

In this paper we discuss real-time stochastic quantization for scalar field 
theory and SU(N) pure gauge theory relevant for QCD. Similar to what 
has been observed in Ref. [U] for scalar fields, also for the gauge theory we 
find that previously reported unstable dynamics represents no problem in 
practice. A combination of sufficiently small Langevin step size and the 
use of a "tilted" real-time contour leads to converging results in general. 
Our procedure respects gauge invariance and appears to be well under 
control. This is exemplified for SU(2) gauge theory in 3 + 1 dimensions. 
For the scalar theory we consider, in particular, a zero-spatial dimension 
example, namely a quantum anharmonic oscillator. The latter is solved 
in addition by alternative methods using Hamiltonian diagonalization for 
comparison. We find that stochastic quantization can accurately describe 
the time evolution for lattices with sufficiently small real-time extent. This 
concerns nonequilibrium and equilibrium simulations at weak as well as 
strong couplings. However, when the real-time extent of the lattice is 
enlarged the stochastic updating does not converge to the correct solution. In 
particular, for the nonabelian gauge theory the thermal equilibrium solution 
is only approached at intermediate Langevin-times. Its life-time in the 
course of the Langevin evolution becomes longer if the complex time path is 
deformed towards Euclidean space-time. 

Since there appears more than one fixed point to which the Langevin flow 
can converge, the solutions obtained by real-time stochastic quantization 
are not unique in general. Similar observations have been made before 
for Euclidean theories with complex actions ^2]- Euclidean theories with 
real actions can be shown to have a unique solution based on positivity 
arguments jl] . A similar argument fails for real-time stochastic quantization. 
In general, here the correct fixed point cannot be chosen a priori without 
implementing further constraints. We prove that all possible fixed point 
solutions fulfil the same infinite set of (symmetrized) Dyson- Schwinger 
identities of the quantum field theory. However, solutions of Dyson-Schwinger 
equations without specifying further constraints are not unique in general. 
We discuss how tests for differentiating them can typically be devised and 
applied. 

The paper is organized as follows. Sec. previews stochastic quantization 
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in Euclidean space-time, which fixes the notation for the derivation of the 
real-time dynamics in Sec. El Particular emphasis is put on the developments 
for gauge theories, since previous suggestions for Langevin dynamics in 
the additional "time" dimension violate Dyson-Schwinger identities of the 
original quantum theory In Sec. 0] we prove Dyson-Schwinger equations 
to follow from real-time stochastic quantization. Simulation results and 
precision tests are discussed in Sees. and El We end with conclusions 
in Sec. 

2 Euclidean stochastic quantization 

The stochastic quantization algorithm calculates ensemble averages in an 
extended space of variables, where the lattice action S of the (^-dimensional 
quantum theory plays the role of a potential energy in a classical Hamiltonian 
H of a (d + l)-dimensional, embedding theory. The ensemble averages of 
the original quantum theory are computed from classical evolution in the 
additional "time" dimension. 

2.1 Euclidean scalar theory 

We first consider classical dynamics for scalar fields in rf-dimensional 
Euclidean space-time with Hamiltonian (see e.g. [13]) 

H= Jd d x Qvr 2 (x; t s ) + C E (^(x; t 5 ), d x ip{x; t 5 ))) . (1) 

Here the field if and momentum ir depend on the <i-dimensional Euclidean 
space-time variable x^ (// = 0, . . . 3) and on an additional "time" variable t§. 
The function C E depends on the field and its derivatives with respect to 
x^. It is taken to correspond to the Lagrange density of the physical four- 
dimensional Euclidean field theory with mass m and quartic self-interaction 
A where 

£e = \ { (^(x; h)) 2 + m\\x- 1 5 ) } + ^ 4 (x; t B ) . (2) 

Since £ contains no field derivatives with respect to t 5 the variable enters 
only as an additional field label in Eq. (|2j). 
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Expectation values of observables F{ip) are given by the functional 
integral 

(F(<p)) = Z~ x J VirVtp F{<p) e~ H[ ^ (3) 
with the normalization 

r VnVtp e - H[ *' v] . (4) 



The Gaussian integrals over the momentum fields can be trivially performed 
and the overall constants cancel out in the computation of field expectation 
values. Here the introduction of the canonical field momenta is used to 
compute expectation values of quantum fields from classical Hamiltonian 
dynamics. This amounts to replacing the canonical ensemble averages (jSJ) 
by micro-canonical ones. The latter can be obtained from solving Hamilton 
equations. 

The basis for the numerical simulations are the discretized classical 
equations for the Hamiltonian (Jl}. The classical dynamics of the field <p and 
its conjugate momentum 7r in t 5 -time is described by Hamilton equations of 
motion 

d(p(x]t B ) 5H 

■ 7r(a?;t 5 ), (5) 



dt 5 5tt(x; t 5 ) 

dn(x;t 5 ) SH 5S E [ 



dt 5 5(p(x;t 5 ) 5(p(x;t 5 



(6) 



where the functional differentiation with respect to ip is for fixed t§. The 
Euclidean action is 

S E [<p] = J d d x£ E ^(x),dMx))- (7) 

In order to discretize Eqs. (jSJ) and © in t$, we denote the difference in 
subsequent i 5 -time steps by At 5 = (t 5 ) n+1 — (t 5 ) n . Correspondingly, we will 
write W((h) n+l ),n((t 5 ) n+1 )} -> {<p', tt'} and W((h) n ),ir((t 5 ) n )} -> {(p,n}. 
A suitable second-order discretization of (J3J) then reads 



= tf*)+<r(,)*.-i^A<i. (8) 
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Correspondingly, expanding n f — n to order At 5 and using Eq. (jBJ) the field 
momenta are (leapfrog) evolved with 



which yields reversible and area preserving discretized dynamics [T3j . 

The conjugate momenta in the classical Hamiltonian (0) have a Gaussian 
distribution, which is independent of the values of the field variables. If 
instead of stepping along a single classical trajectory the momenta after every 
single step are randomly refreshed, one recovers classical Langevin dynamics 
in £ 5 -time. It amounts to making the substitutions 

V2tt(x) - -n{x) , ~A<*-e, (io) 

with Gaussian noise 



(^(ar)), = , (rj(x)rj(x')) T] = 2 5(x-x'), (11) 
where the average of an observable A((p; rj) over the noise is given by 

_ J[di)M(y;i))exp{-l S& l xrf(x)} 

{ ,)> " = jw^Fw^m ■ (12) 

The discretized Langevin equation then reads 1 

tf(x) = <p(x) - + y/ir,(x) . (13) 

The sum over the Langevin steps, d = ^2 e, is called " Langevin-time" , which 
replaces t 5 of the classical dynamics. 

This simplified description of the classical dynamics forms the basis of 
stochastic quantization The stochastic process (|T3*j) is associated to 
a distribution Pe(^s) for the field ip(x). Writing Pe( (£5)71+1) — > Pr an d 
-Pe((^5)ti) — ► Pe its evolution can be obtained from 

PkW\ = ( / [d^PMllS (y(x) - <p(x) + - yfa(x)j } . (14) 

1 Our discretization corresponds to Ito calculus for stochastic processes. 
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Expanding the 5-functionals and keeping only terms up to order e gives the 
Fokker-Planck equation 

= /*d««A(^ + fl ,«|) M + 0(e) , (15) 

where we have used Eq. (fTTj) . For real action S E one can prove that it 
converges to the late Langevin-time limit 

Pe[<p] =e~ s ^ + 0(e). (16) 

Therefore, expectation values (F((p)) can be obtained from noise averages 
or, assuming ergodicity, from Langevin-time averages for sufficiently long 
classical trajectories. 



2.2 Euclidean gauge theory 

A similar discussion can be done for gauge theories. We consider a nonabelian 
pure gauge theory on an anisotropic lattice of size (N s a s ) 3 x N T a T with 
Euclidean action 

Se[U] = -^EEEl^T^^ + T^i)- 1 } 

( Tr u ^ + Tr u £) - 1 } > ( ir ) 



i<3 



with spatial indices i,j = 1,2,3. It is described in terms of the gauge 
invariant plaquette variable 

U x ,y,u = ^ x,pP x+p,,vU x+i/,^ x,v , (18) 

where U~l = U X;IXU . Here U XtfJ , is the parallel transporter associated with the 
link from the neighboring lattice point x + p, to the point x in the direction 
of the lattice axis /i = 0, 1, 2, 3 with = _^ • The link variable U x>fi is 
an element of the gauge group G. Because of the anisotropic lattice we have 
introduced the anisotropic bare couplings g for the time-like plaquettes and 
g s for the space-like plaquettes with 

(19) 

9o 9,1e 
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where 7e = a s /a T is the anisotropy parameter. 

For G = SU(N) one has U~ x = U]. . When we consider Minkowski 
space-time below we will observe that the latter no longer holds. Therefore, 
we keep U x \ v in the definition of the action (JT7J), which will still be valid for 
the Minkowskian theory. 

In order to derive the ts-time dynamics of the 5-dimensional theory one 
can follow the equivalent steps as in Sec. 12.11 For the definition of the 
Langevin drift term one has to define differentiation with respect to the 
nonabelian variable U Xjll . Differentiation in group space will be defined by 

D xlxa f{U X)il ) = A/ (e iaX «U x ^) U =0 (20) 

with the generators A a of the Lie algebra. For G = SU(N) these are the 
a = 1, . . . , N 2 — 1 traceless, hermitian N x N Gell-Mann matrices, which are 
normalized to Tr(A a A&) = 25 ab with [A a , A&] = 2if abc \ c , where the structure 
constants f a b c are completely antisymmetric and real. The derivatives then 
satisfy the commutation relations [D Xfia , D Xflb ] = 2f abc D xfJ/C . 
For the Hamiltonian 

H = \Y, P la{h) + S^U{t,)} (21) 

xfia 

with the above definitions the Hamilton equations for the ts-time dynamics 
read 

o, ' = ^aPx/iaUx^ , —7- = — D xf j, a £>E[U] ■ (22) 

The discretized canonical equations of motion corresponding to (JHJ) and Q 
are then given by 



dU XjlJb . ld 2 U x 

^r At5+ 2^i 



1 d 



■ r/iu 1 ' j-./i ;> T 2 Q£ \ y o, Xfia'-J x,fj,J '-» t '5 



U x ,n ~\~ i^aPx/j,aUx,n A^5 + n 0j (yi^a-Px(j,aU x ,i_i) At 

exp <j i\ a [ P XIM At s - ^D XIM S E [U\ At 2 ) } U Xttt , (23) 

where the equalities hold up to corrections 0(At|). Here P" is the t 5 -time 
conjugate momentum for U x ^ with (leapfrog) discretized Hamilton equation 
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(cf. ©) 

Km = P x , a - I {D Xfia S E [U] + D' xlia S E [U'}) At 5 , (24) 

where D' refers to the derivative with respect to U' . 

The Langevin equation is obtained in the same way as described for the 
scalar case by the corresponding substitutions (jTUjl with 

V2P X 

fia 

(x) — > T) XIMl (x) , (25) 
and a real Gaussian noise rj Xfia satisfying 

{Vxfj,a)r] = 0, (j)xiia.T)yvb) T) = 2 Sf ll/ S X yS a b . (26) 

The discretized Langevin equations for U Xjlx to this order in e may then be 
written as 

u L,n = ex P {-^a (eD Xfia S E [U) - y/er} xlu ) } U Xt(t . (27) 

3 Real-time dynamics 
3.1 Real-time scalar theory 

Instead of the embedded <i-dimensional Euclidean theory discussed above, 
one may consider Minkowskian space-time. This requires an analytic 
continuation of the Euclidean time in the action (JJj) to Minkowski time. 
For this it is instructive to consider the one-parameter family of actions S% 
for the field <p{x) with x = (t, x) given by 

S&] = ~ J d d xe-^ 2 j^D^ + imV 2 + , (28) 

with d'Alembertian 

□ { = e l ^d 2 t - V 2 . (29) 

For £ = one recovers from Eq. ()28j) the standard Minkowski action, while 
for £ = 1 one has 

S^=i = i S E • (30) 
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The discretized Langevin equation (J13|) can be written for the family of 
actions (J28j) as 

ip '( x ) = ip ( x ) + iep^- + V~ev(x). (31) 
bip{x) 



with 

^ = {rk<p{x) + m*<p{x) + ^) ) • (32) 

The Gaussian noise r](x) is defined as in the Euclidean case in Eq. (fTT|) . The 
case £ = forms the basis of real-time stochastic quantization [Hj, and we 
will denote S^ =Q = S. 

It is important to note that possible solutions of Eq. (j3*Tj) will not be 
real in general. For instance, for a real scalar field theory with Minkowskian 
action S[(p] it will generate complex field values for > 0. For a complex 
field 

<p(t, x) = (p R (t, x) +i(fi(t, x) (33) 
also the conjugate momenta 7i(t, x) or, with (JTUJ) , the respective noise 

r](t, x) = r] R (t, x) + ir]i(t, x) (34) 

can be complex. Equation (pTT]) may then be written as 

</4(^ x ) = *) - eI^(p R ,ipj\t,x) + y/e7} R (t,x) , 

x) = *) +eR^(cp R ,(p I ;t,x) + y/er]i(t,x) , (35) 

where 

i? c (^,^;t,x) = Re( /^ ] ) , 

\dip\t, Xj tp=<pa+i>pi J 

I^, W *,x) = Imf-^L . (36) 

\ <)<£>(£, X) tp=<PR+itpiJ 

It remains to consider suitable choices of the noise terms. A particularly 
simple choice is a real noise, where rjj = 0. However, for complex fields a 
real noise may not be suitable in general, since the noise plays the role of the 
field conjugate momentum in the underlying classical dynamics (see (jlOjO . 
For a complex noise one obtains from (|11|) 2 



2 Note that the complex field ip yields a non-real action S[tp], which is not a functional of 
if* p. Likewise, the conjugate momentum or complex noise does not fulfil (r)*(x) rj{x')) v = 
28{x-x') but HU). 
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(Vr(x) rj R (x')) v = 2a5(x - x') , (rji(x) r] I {x')) v = 2(a - l)6(x - x') . (37) 
For real a > 1 the average of an observable A(<p; rf) over the noise is defined 



as 



/[d77 fl ][d77j] Aexp {- Jd 4 x (-L t£ + rf) } 

(A) v = i ^ y J • (38) 

/[d^fd^] exp |- Jd*x r/l + _L_ | 

The case a = 1 corresponds to a real noise with rji ee and noise average as 
in (I12|) with r\ — > 77^. We will see in Sec. EJthat the different possible choices 
do not affect the final result. 

Other observables F((p(x)) follow similar evolution equations, which are 
obtained like (jHJ) from discretization to second order in At§, 



r,, // xx „, / u dFUpix)) . 19 2 F(oj(x)) . 2 
F ifpfix)) = F ^(x)) + ^ Ats + 2 A ^ ' (39) 



with the substitutions (jTOj) and using (|31|). E.g. the product of two fields 
follows the discretized equation 

<f/(x)<f/(y) = ^x^ + ie^^^ + ^M^x) 

+ erj(x)r](y) + \Te {if(x)r](y) + ip(y)r](x)} . (40) 
For three powers of the field one finds 

<p'(xW(y)<p'(z) = <p(x)<p(y)<p{z) 

+ y/e {<f(x)ip(y)r](z) + if(x)ip(z)r](y) + <p(y)ip(z)r](x)} 

- e { fM v{yMz) + s m v(xMz) + fM v(xMy) 

I 0(p[x) btp{y) dtp{z) 



e 



{(p(x)rj(y)r](z) + (p(y)ri(x)ij(z) + <^{z)r](x)r](y)} , (41) 



and correspondingly for higher powers of the field. In Sec. H]we will see that 
the late Langevin-time limit of these equations leads to an infinite hierarchy 
of Dyson-Schwinger identities. 

3.2 Real-time gauge theory 

Following Sec. 13. II we replace the Euclidean action (fTTj) by the Minkowskian 
with 

- S E [U] = iS[U) (42) 
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on a real-time lattice of size (N s a s ) 3 x N t a t . The real-time classical action 
reads 

S[U] = -^^^/-^-(Tr^ + TrCZ-y-lj 

+AEE{m ( Tr ^ + Tr ^i) - 



i<3 



(43) 



where the relative sign between the time-like and the space-like plaquette 
terms reflects the Minkowski metric, and 

As *z£i, As ™, (44) 

with the anisotropy parameter 7 = a s /a t . 3 

The discretized Langevin equations for U Xtft to second order in e may then 
be written as 

U' X41 = exp {i\ a (e iD a . A40 S , [E/] + y/erj Xfia ) } U XtlM , (45) 

with 

iD xtia S[U] = T^yPo ^2 fe Tr {^aU x jC xJ0 ~ C xJ0 U~jX a ) 

j 

+^0AiTr (\ a U X fiC X fij — C X fijU X Q\ a ) } 



2N 

hi 



I 3 

- /V* (*aU x ,»C x ^ - C x ^U x ,lXa) ■ (46) 



2N 



For a compact notation we have defined (5^ = (3 S and (3 0i = /3 i0 = — A) an d 
- TT ~ TT~ l TT~ l 4- TT~ l TT~ l TT - 

C x ^ u — U XjV U x +& } ijU x _^^ v + U x _~^ I/ U x ^£,^U x +f l ^0 il/ . (47) 



5 See Sect 12. 21 for a t — —ia T and the replacement according to Eq. Q42J1 . 



11 



With U X ^C X ^ U = U Xtia ,+U Xi n(- v ) and C Xtltu U x * = U^+U^^ one observes 
that the sum in Eq. (pftjj) is over all possible plaquettes containing U XjfJi . 

Specifying to SU(2) gauge theory A a = a a (a = 1,2,3) represent the 
Pauli matrices, and we can make further simplifications using Tr([/ _1 cx a ) = 
— Tr(Ua a ) for any element U G SU(2). The latter simplification also holds 
for U G SL(2, C). This is relevant since, similar to what has been discussed 
for the scalar theory above, possible solutions of Eq. may respect an 
enlarged symmetry group. Taking 

U Xill = e iAx » aaa/2 (48) 

the vector fields A xtia need not to be real, which is in contrast to the Euclidean 
case discussed in Sec. 12.21 The complex matrix A a xl) a a still remains traceless, 
however, the Hermiticity properties are lost. As a consequence, it is no longer 
possible to identify W with £7 _1 as is taken into account in Eq. ()17|). This 
corresponds to an extension of the original SU{2) manifold to SL(2, C) for 
the Langevin dynamics. Only after taking noise averages the original 577(2) 
gauge theory is to be recovered (see Sec. 0J). 
According to Eq. (j2lj|) the noise is given by 

{Vxfia)rj , (Vxfia ^]yvb)r\ 2b^ lv b X yb a \ ) . (^9) 

It is essential to use the same statistics for the noise as one has in the 
Euclidean case. One may be tempted to replace the b^ u on the right-hand side 
of Eq. by g^ v to make the Langevin equation manifestly covariant 0. 
However, in this case solutions of the Langevin evolution would not respect 
the Dyson-Schwinger identities of the underlying quantum field theory, as 
is shown in section Only for observables one has to respect Lorentz 
symmetries which is still fulfilled with Eq. (149)) . 

In general, observables F(U) follow similar evolution equations, which 
are obtained from the second-order discretization equivalent to Eq. ()39|). For 
instance, for the plaquette U X)IU/ given by (fTHj) one has 

= u ^ + ^ + Atl ' (50) 

with Eqs. ()22|) and ()4*2*|) . This will be used to derive the corresponding Dyson- 
Schwinger equation for the plaquette variable in Sec. 14.21 
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4 Fixed points of the Langevin flow in real 
time 



4.1 Scalar theory 

Starting from some (arbitrary) initial field value <f(x) at tc, = we compute 
the Langevin-time flow according to equation (J3TJ. A fixed point of 
the evolution equation is defined by the condition for the noise-averaged 
field ((p'(x)) v — ((p(x)) v = 0. It simultaneously corresponds to a fixed 
point in the space of all noise-averaged correlation functions (^^^(y))^, 
(ip(x)ip(y)ip(z)) ri , etc. following equations (J1U|) . PT|) and similarly for higher 
n-point functions. According to Eq. (J31j) a fixed point for the one-point 
function corresponds to 



/5S< 



W( 



X 



0, (51) 



where we have used Eq. (II lj) or equivalently Eq. (|37|). It is important to note 
that the result is independent of the different possible implementations of the 
noise according to Eq. (}3*Tj) . The two-point function described by Eq. (jlUj) 
then fulfils 

(wH^lMr 2 *^' (52) 

For the 3-point function one has from (|41l) 



/ r? 



2i ({y{x)) n 5{y - z) + {y{y)) v 5{x - z) + ( V (z)) v 5{x - y)) , (53) 

and correspondingly for the higher n-point functions. Equations (|5*T|) - (|5!3|) 
and the corresponding equations for n > 3 are the symmetrized Dyson- 
Schwinger identities for the time-ordered correlation functions of the scalar 
theory in Minkowski space-time. 
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4.2 Gauge theory 

For the plaquette variable (|5U|l one has 

d 2 U x , t iv \ / d 2 U x ^ x x d U x +p,,v TI -x jj-i 

J \ Q£2 u x+fi,u^ x +v,ti U x,v ' u x,[i U x+0, l i U x,v 

-I- 77 77 - X+V^ TT-l ,JT JT _ 7-T- 1 

5 5 

/ n 

= 0, (54) 

where we have used that noise averages over products of first-order 
^-derivatives vanish. The latter can be observed from the fact that link 
and momentum variables are not correlated such that their noise average 
factorizes. According to Eq. ()49|) the averages over the corresponding 
momenta vanish identically. The symmetrized Dyson-Schwinger equation 
for the plaquette variable then follows from Eqs. (J22J), (J4~2|) and ()25|). i.e. 



d 2 U x ^ 

—7^ = 1 



'5 

d 2 U x ,l 
dt 2 



. )f2 - Y ^ a i iD ^aS[U}) U X4L - i Y KVxfjahVx,jbU Xill , 



a a.b 



^bVx^b , (55) 



with Eq. and 

Y ( Aa )«/3 ( Xa \s = 2 (tasSpj - j^6 a p8 lS j • (56) 

By setting the first term of the sum in Eq. (|H%j) to zero and taking the trace, 
one obtains, using the notation of Eq. (jlRjl : 



jy 2 i ^ / 

2 — (TrL^r^j/)^ = — ^ ^ /3^7 \ Tr ( y C x ^~ 1 U x ^U x ^ l , U x ^C x ^-yU X:fMU j 



7=0 
7#M 



— ^Tr£4^Tr (C Xtfll U x j l — U X ^C X ^) \ (57) 

The first term on the RHS of this equation contains U X)fi and its inverse 
such that the loop can be viewed as including a departure in the 7-direction. 
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Figure 1: Graphical representation of the Dyson-Schwinger equation for a 
Wilson loop. 

The second term contains U Xtft twice: The loop can be viewed as including 
a "crossed path" (see Fig. The last terms contain separately traced 
plaquettes. For U G SL(2, C) or SU{2) one has TrU = Tr£/ -1 such that 
these last terms in Eq. (|F7j) vanishes, i.e. Tr {U X ^C X41U — C x> ^U~t) = in 
this case. A similar analysis can also be done for the second to fourth terms 
of Eq. ([54)1 . such that the sum is symmetric in all the links of the plaquette 
U X)fll/ . In Fig.[T]the Dyson-Schwinger equation is displayed graphically, where 
it is generalized to an arbitrary closed loop. When a certain link appears more 
than once in the original loop then additional "contact" terms appear. The 
respective analysis follows along the same lines and will not be considered 
here. 

Fixed points of the Langevin flow for observables F[U] representing gauge- 
invariant products of link variables along closed loops are obtained from 
(d 2 F[U)/dtf) v = 0. The corresponding set of Dyson-Schwinger identities 
[T4*j constitute an infinite system of equations whose solution is not unique 
in general without specifying further boundary conditions. Accordingly, we 
will typically find more than one possible fixed point below, when we solve 
the Langevin equation numerically. They all solve the same set of Dyson- 
Schwinger equations. In contrast to the case of Euclidean space-time, where 
the associated Fokker-Planck equation (|THj) can be shown to converge to a 
unique solution at late Langevin-time, there seems no general proof for real 
times, i.e. in the absence of a positive definite probability weight. 
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5 Accuracy tests I: scalar theory 



In thermal equilibrium the fields obey the periodicity condition ip(0, x) = 
(p{— z/3,x) with inverse temperature (3. Accordingly, correlation functions 
in Euclidean space-time can be computed using a purely imaginary time- 
path from t — to — if3. Thermal equilibrium correlation functions 
(ip{xi)ip(x2) ■ ■ ■ (fixn)) with real times £i, £2, • ■ • t n have to be computed using 
a time-path that extends along the real-time axis including these times. 
The curves on the left of Fig. |21 give some examples of possible real-time 
contours for thermal equilibrium along with other complex contours that will 
be employed below. The curves (also denoted as "closed" and "rectangular") 
both first run along the real-time axis and then turn in different ways to 
—i(3. The curves "isosceles" and "asymmetric" exhibit a tilt with respect to 
the real-time axis. For nonequilibrium evolutions we will use the isosceles 
triangle or asymmetric time-path below, where a small tilt with respect to 
the real-time axis can serve as a regulator to improve convergence. 

We consider field theory for the different choices of complex contours 
shown in Fig. |21 We discretize the contours using the complex-time points Ct 
with t — 0, . . . , Nt and ReCo = Re(7/v t . With the notation <£>t(x) = <p(Cf, x) 
and A t = C t +i — C t a scalar theory on such a contour is defined by the action 



S = - > ' / d 3 x 



1 



(V?t+i(x) - ¥?i(x)) s 



2 ^ J 1 At 



+Y [v?m(x)VVm(x) + ^(x)VVt(x) 
-A t [V(<p t+1 (x)) + V(<p t (x))}\. (58) 



It is instructive to consider for a moment the free theory with V((p) = m 2 ip 2 /2 
and neglect spatial dimensions for simplicity. For the free theory the action 
(|58|) is quadratic in the fields: 

s = iX> G -V'< • = (59) 

The complex inverse contour propagator G^} is symmetric and not Hermitian 
in general. With 

J2G;M = c a ^ (so) 
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Re 



Figure 2: Complex time-contours (left) and the corresponding distribution 
of eigenvalues (right) characterizing the Langevin dynamics for a free scalar 
theory as described in the main text. 

its eigenvectors iff and complex eigenvalues c a can be used to write the action 

as 

S = \Y. ca ^ a X a , X° = I>iVt> (61) 

a t 

where we have used completeness and orthogonality relations Et iff iff — $ ab 
and *£ a if?if?, = tw- 
in Fig. |21 we show the distribution of complex eigenvalues c a for the 
different corresponding contours displayed on the left. The eigenvalues can 
become small though they are non-zero. They depend strongly on the chosen 
contour. Its impact on the Langevin dynamics can be observed from equation 
(pTT|). which becomes here a set of independent Langevin equations 

x af = x a + iec a X a + V^V a , (v a V b ) = ^ ab , (62) 

with x a — Et i^tVt- F° r a purely imaginary time-contour from zero to —if} 
one has Rec a = 0, Imc a > (Va), as for the Euclidean case described by 
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Eq. The Langevin evolution in £ 5 -time converges for Imc a > jH] and 

one finds in the continuum limit 



(xY>„ + 



26' 



(63) 



i(c a + c b ) 



Therefore, eigenvalues with small positive imaginary part can converge 
slowly. However, this linear or free field theory analysis becomes invalid in 
the presence of sufficiently strong interactions. The associated damping due 
to interactions in realistic field theories typically leads to rapid convergence, 
as can be observed from the numerical results below. (See also Ref. [Hj.) 

Fluctuations can grow large if the dynamics is governed by small real 
and imaginary part of the eigenvalues c a . The linear analysis indicates that 
this should be a problem for contours along the real-time axis. From the 
eigenvalue distribution of Fig. El one observes that a non- vanishing tilt of the 
contour with respect to the real-time axis may serve as a regulator. For our 
simulations we employ a small tilt such that the contour always proceeds 
downwards, i.e. \m.C t > ImC t+ iVt. Without or for very small regulator 
we encounter incidences of unstable Langevin dynamics (see also [12J ) . 
Their appearance depends on the random numbers and they are strongly 
suppressed by using a smaller Langevin step size. We normally discard these 
trajectories and typically employ a Langevin step size of e ~ 1CT 5 . 

5.1 Short-time evolution: Thermal fixed point 

We consider an interacting scalar theory with classical action (JSHJ), where 



in zero spatial dimension. The Schrodinger equation for the corresponding 
quantum anharmonic oscillator can be solved numerically by diagonalization 
of the Hamiltonian. We use the occupation number representation and 
truncate the Hilbert space keeping 32 dimensions corresponding to the lowest 
occupation numbers. Exponentiation of the diagonalized Hamiltonian yields 
the time translation operator, which may be used for real as well as complex 
times. We checked that the results are insensitive to an enlargement of 
the truncation dimension. This is compared to the results obtained from 
stochastic quantization by solving the Langevin dynamics according to (|31|) . 



V(<p) = \m\* + ±<p* , 



(64) 
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Figure 3: Real and imaginary part of the unequal-time two-point correlation 
function for the quantum anharmonic oscillator as a function of real time 
t. The results obtained from stochastic quantization agree to very good 
accuracy to those obtained from directly solving the Schrodinger equation. 

In the following all quantities will be given in appropriate units of the 
mass parameter m. We first perform simulations in thermal equilibrium, in 
which case the time contour extends along the imaginary axis to the inverse 
temperature (3 (see Fig- EJ) - Non-equilibrium is considered in Sec. I5.3I 

Fig- Elshows the two-point correlation function (</?(0)</?(t)) T? as a function 
of real time t for A = 24 and (3=1. The real-time extent of the contour 
is tfinai = 0.8, such that the upper branch has a tilt of 0.01/9 and the lower 
branch has a tilt of 0.99/5. Square symbols denote the results from stochastic 
quantization for the real part of the two-pint function and crosses for the 
imaginary part, respectively. For comparison the corresponding results from 
the Schrodinger equation are given, which agree well. 4 

Fig. |3] shows the real part of the two-point correlation function as a 
function of real time t, where we employ a much stronger coupling A = 96 
and (3 = 1. Square symbols denote the results from stochastic quantization. 

4 For Figs.GHElthe total number of points along the contour ranges from 32 to 64, with 
equal number of points on each branch. The noise-average is obtained from the average 
over 10 3 to about 10 4 runs. The Langevin step size is e ~ 10~ 4 to 10~ 5 . 
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Figure 4: Similar evolution as in Figs. El but for much stronger coupling 
A = 96 and (3 = 1. 

Here the real-time extent of the time contour is ifi na j = 0.5. The upper 
branch of the contour has a tilt of 0.001/3 with respect to the real-time axis, 
so it is almost horizontal and, therefore, realizes to high accuracy a real- 
time contour. For comparison results from the Schrodinger equation are 
displayed as well. The solid line corresponds to Minkowski time evolution 
("real contour"), whereas the dashed line gives the results for those complex 
times with small imaginary part as employed in the stochastic quantization 
simulations ("complex contour"). They all agree well. We fit the time 
evolution to 

(<p(0)<p(t)) = a cos (w Ret) e" 7Re * , (65) 

which allows us to extract characteristic damping times. We find from the 
stochastic quantization simulation = 0.018 and from the Schrodinger 
equation j/u — 0.013 (Minkowski) and 7/0; = 0.014 on the contour. The 
frequency is always uj = 2.95. 

Fig. shows a time evolution for smaller temperature with j3 — 8 and 
A = 6. The real-time extent of the contour is tg na i = 4. Using the fit (J55J) 
we find from the stochastic quantization simulation 7/0? = 0.92 and from 
the Hamiltonian diagonalization method 7/cj = 1.0 on the contour. The 
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Figure 5: Similar evolution as in Figs. El but for smaller temperature ((3 = 8 
and A = 6) on an isosceles triangle contour. 

frequency is u = 1.4. 

Fig. |U1 verifies the validity of the Dyson-Schwinger identity ()52j) for various 
real-time values. We use the parameters A = 24 and (3=1. The real-time 
extent of the contour is tfi na i = 1 with 32 points in total on the (isosceles 
triangle) contour. Plotted is the Langevin-time evolution of the LHS and 
the RHS of the equation 5 

Yl G oh^wt>) - 8tt> = -i^(<pm<pm>) , (66) 

t 

where Gq * f is the inverse propagator (}59|) of the free theory. One observes that 
at sufficiently late Langevin times both sides agree to very good accuracy. 
In particular, one sees that equal-time values at different real times agree, 
which has to be the case for the time translation invariant thermal solution. 



5 Eq. (|52() is the symmetrized form of (|66|l . 
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Figure 6: Shown is the Langevin-time evolution of the two-point function 
and four-point function given by the LHS and RHS of the Dyson- Schwinger 
identity (jOfcij) . The curves are for fixed real-time values. At sufficiently late 
Langevin-time they agree to very good accuracy, thus becoming related by 
the Dyson-Schwinger equation. 

5.2 Long-time evolution: Non-unitary fixed points 

We observed above that stochastic quantization can describe the real-time 
evolution very accurately for short times. However, when the real-time extent 
of the lattice is enlarged the stochastic updating may not converge to the 
correct solution. As an example, Fig. [7| shows the two-point correlation 
function as a function of real time t for A = 24 and (3=1, i.e. for the 
parameters of Fig- El From that figure we have seen that for a real-time lattice 
with tfinai = 0.8 there is excellent agreement of stochastic quantization and 
Schrodinger equation results. Correspondingly, Fig. [7| exhibits a constant 
real part and a vanishing imaginary part of the equal-time correlator in 
thermal equilibrium for t final = 0.8. However, doubling the extent of the 
contour leads to a qualitatively different behavior as indicated by the cross 
symbols of Fig. EJ This difference persists also on finer real-time grids. The 
non-vanishing imaginary part of the equal-time correlator and the loss of 
time-translation invariance reflects a non-unitary time evolution. 
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Figure 7: Real and imaginary part of the equal-time correlator as a function 
of time (index of lattice site along the time-contour) for the parameters 
(A = 24 and (3 = 1) as well as the contour of Fig. |3J Compared are two 
simulations where the real-time extent of the lattice differs by a factor of two. 
The larger lattice leads to a qualitatively different, non-unitary behavior. 

The properties of these non-unitary fixed point solutions depend on the 
details of the time contour. This is in contrast to the universal properties of 
the thermal solution. As an example, Fig. |H] shows the equal-time two point 
function for same A = 24 and (3 — 1 as in Fig. [7[ however, with different 
contour geometry (isosceles triangle). For t&aa\ = 1 one still observes the 
thermal fixed point solution, whose properties agree very well with those 
obtained from lattices with smaller £fi na i- Stretching the temporal extent of 
the contour to tfi na i = 2 the Langevin dynamics converges to a non-unitary 
fixed point as one can see from Fig. |H| The detailed properties of this fixed 
point differ from those shown in Fig. [7| 

In Fig. |H1 we have shown how the propagator fulfils for the thermal fixed 
point the corresponding Dyson- Schwinger identity. According to the results 
of Sec. H]any fixed point solution has to fulfil the same properly symmetrized 
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Figure 8: Similar evolution as in Fig. [7[ but on an isosceles triangle contour. 
While the thermal solution obtained for t final = 1 does not depend on the 
contour details (see also Fig. [HJ), the properties of the non-unitary fixed point 
obtained for i final = 2 are contour dependent (see also Figs. El and [TUJ). 

Dyson-Schwinger identities. For instance, the identity reads 
it 

(67) 

We display the LHS and the RHS of this equation as a function of Langevin- 
time in Fig. Indeed, one observes that the symmetrized Dyson-Schwinger 
equation (|52j) holds within statistical errors at late Langevin-time for the non- 
unitary solution. In Fig. ^] we check Eq. ()66|). which is the non-symmetrized 
version of Eq. (|52|). Here the discrepancy we find for the non-unitary fixed 
point is greater than the statistical error. In contrast, the thermal fixed 
point solution respects, of course, both Eqs. (jBTjl and as required by any 
physical solution for the underlying quantum field theory. 
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Figure 9: The left-hand-side (LHS) and right-hand-side (RHS) of the Dyson- 
Schwinger equation (JBTj) for the non-unitary fixed point, which agree well at 
late Langevin time. Contrary to the properties of the thermal fixed point 
(see Fig. EI), the (0,0) and (t,t) components are not equal, which reflects the 
loss of time translation invariance. 

5.3 Non-equilibrium dynamics 

Non-equilibrium dynamics can be described by the generating functional for 
correlation functions [To] : 

Z[J;p] =Tr{pT c e^(-)*Mj 
= J d<pid(p 2 p((pi,(p2) J Z)Ve !/c(iW+J(xMx,) . (68) 

The path integral (jUS|) displays the quantum fluctuations for a theory with 
Lagrangian L, and the statistical fluctuations encoded in the weighted 
average with the initial-time (non-thermal) density matrix p(<^i, y^)- Here 
Tc denotes contour time ordering along a closed path C starting at t = x° = 
with j c = j c dx° J d d x. This corresponds to usual time ordering along the 
forward piece C + , and anti-temporal ordering on the backward piece C_. 
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Figure 10: The left-hand-side (LHS) and right- hand- side (RHS) of equation 
(HSU), showing a discrepancy beyond the statistical error. 

Denoting fields on C + by tp+(x) and on C_ by y?_(x), the initial fields are 
fixed by fi{x) = <f+(t = 0,x) and y^fx) = <f-(t = 0,x). Non-equilibrium 
correlation functions are obtained by functional differentiation and setting 
J = 0. 

With this notation the expectation value of a real-time observable A(ip) 
can be written as 

J ip + (0)=ipi 

D V _D V+ e iS ^+' lp - ] A{^ + ) . (70) 

Here S p contains the actions on both contour branches as well as the density 
operator. In the following we consider Gaussian initial density matrices for 
which, neglecting all space dependence for simplicity, the most general S p 
reads [12] 

S p [cp + ,cp_] = S[ip + }-S[ip-]--So(<p + (t = 0),cp-{t = 0)}, (71) 

at 

S [<f+,<f-] = i0(v?+ - (p~) - - g ~^ 2 1 ((<p+ - 4>f + (<P- - 4>) 2 ) 
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+ §((,+-0) 2 -(,--0) 2 ) 

+ ^r(<p+-<P)(<p--<P)- (72) 

The real parameters <p, <f), a, £, rj determine a complete set of independent 
initial one-point and two-point correlation functions: 



<i> = (<p(t = o)) , <p = m = o)) 

e = (<p{t=oMt=Q)) e , 
i 

2 



^ = -(^ = oM* = o) + ^ = oMf = o)) c) 



a 2 



^ +4^2 = = Wit = 0))c • (73) 

Starting from a given initial density matrix the nonequilibrium simulation 
is carried out using the Langevin equation (|3T|) with the action replaced by 
(HQ), i-e. 

V4(x) = y» ± (x) + i 6 + v^^(x) , (74) 

bip±(x) 

updating all the points including <p±(t = 0). 

As an example, Fig. El shows the time evolution of the expectation 
value (<f(t)) as a function of time, all in units of the mass parameter m 
of equation ([640. Here the coupling is A = 1 for zero spatial dimension. 
The initial conditions are (ip(t = 0)) = 1, (ip(t = Q)ip(t = 0)) c = 1, 
(cp(t = 0)<f(t = 0)) c = 0.25 and zero for the remaining quantities of (J73J). 
The time contour is tilted with 5% slope, which is denoted as complex 
contour in Fig. El Results using stochastic quantization are given for three 
different lengths of the real-time extent of the contour (tfi na i = 2, 1 and 0.5). 
For comparison we also show the results of the numerical solution of the 
corresponding Schrodinger equation (see Sec. 15.1(1 for real times as well as 
for the complex times on the tilted contour. This also illustrates the contour 
dependence of the results. One observes that the stochastic quantization 
algorithm produces accurate results on the contour for sufficiently short real- 
time extent. Like for the thermal equilibrium results of Sec. 15. 1( we find that 
it becomes inaccurate for late-time evolution. 
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Figure 11: Non-equilibrium time evolution of the field expectation value 
as a function of time. Shown are results for different real-time extents of 
the lattice. Shorter i final lead to improved agreement with results of the 
Schrodinger equation on the employed complex contour. 

6 Accuracy tests II: nonabelian gauge theory 

In the following we perform a similar investigation for SU(2) gauge theory as 
done above for the scalar theory. Though some aspects of the application of 
real-time stochastic quantization are comparable, we will find that there are 
crucial additional restrictions concerning its range of validity for nonabelian 
gauge theories. 

The equation for the Langevin evolution is given by Eq. PB^l for the 
action (H2J), where we employ the approximation that g = g s . We use the 
U = al+ib a a a , a 2 +b 2 a = 1 representation for SU(2) matrices. For SL(2), the 
representation with the same constraint can be used, but the parameters a 
and b a are no longer real. We expand the exponential to first order in e, which 
means we must include the square of the noise term, which is proportional 
to unity. The evolution equation then reads 

U' x>fM = {al + ia a (eiD x ^S[U] + ^Vx^))U x ^. (75) 

To stay in group space the constant a is calculated from the constraint 
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Figure 12: The spatial plaquette average as a function of Langevin time. 
Shown are results for different complex contours. Here a = corresponds 
to a contour with an infinite extent along the real-time axis, while a = it/ 2 
denotes the Euclidean contour. The longer the real-time component of the 
contour the less accurately the thermal solution is reproduced. 

a 2 + b\ = 1 of the matrix transforming U X)tt . 

As the starting configuration for the Langevin-time evolution we take all 
link variables equal to unity. Typically we employ I/7 = a t /a s = 0.01 — 0.5. 
We use lattices of spatial size N 3 = 4 3 . The Langevin step used is ~ 10 -6 . 
All quantities are given in units of a s . The triangle contour of Fig- EI is used, 
with r + being the imaginary extent of the contour on the forward part and 
t_ on the backward part. We calculate thermal distributions, with inverse 
temperature r + + t_ = 4. 

Fig. [T2] shows the Langevin-time evolution of the spatial plaquette 
average. 6 The solid line shows the result for vanishing real-time extent, i.e. 
for a Euclidean contour. The different dashed curves correspond to results for 
complex contours on isosceles triangles (r + = t_ = 2) each having a different 
tilt a with respect to the real-time axis. Here a = would correspond to 
an infinite extent along the real axis. One observes that with increased real- 

6 The plaquettes are also averaged in real (complex) time. 
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Figure 13: The Langevin-time "# crc .ssover for the crossover away from the 
approximate thermal solution to another fixed point. (See Fig. fT21 ) It is 
displayed as a function of tana for various time contours. Contours with 
different real-time extent, but with the same r+ are connected with a line. 

time extent or smaller tilt a the correct thermal solution is approached less 
accurately. 

In particular, it is only approached at intermediate Langevin-times, 
irrespective of the details of the non-vanishing real-part of the contour. This 
aspect differs from the scalar case where short real-time extents lead to stable 
thermal solutions. For the nonabelian gauge theory the correct thermal fixed 
point is approached at first. However, it is not stable and the Langevin flow 
exhibits a crossover to another (stable) fixed point. The fact that the thermal 
solution is not a true fixed point for non-Euclidean contours can be seen by 
monitoring for example (lm(iTr(U X)fl a a )) 2 ) , which is zero for SU(2) matrices 
in the Euclidean theory. It is found to grow exponentially while observables 
of the original SU (2) theory such as the plaquette average still approach the 
thermal solution. 

This is quantified in Fig. where the L angevin-time ^crossover? ^f about 
which the crossover from the approximate thermal solution to another fixed 
point occurs, is given for various time contours. One observes that the time 
of the crossover is mostly dependent on the angle of the slope of the contour. 
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Figure 14: The spatial plaquette averages as a function of the index of 
lattice sites along the time-contour. One observes that it is time translation 
invariant to rather good accuracy for the approximate thermal fixed point, 
while it is not for the non- unitary fixed point. 

Connecting points with the same r+, one sees that Crossover is approximately 
proportional to the logarithm of the tangent of the tilt of the contour. Here 
"^crossover was measured on N t = 16 lattices, using 5 runs per parameters. 

In Fig.Elthe spatial plaquette average is shown as a function of the index 
of lattice sites along the time-contour. While the Langevin-time evolution 
approaches the thermal solution it is time translation invariant to rather good 
accuracy. The observed small discrepancy is decreasing with Langevin-time 
while the system stays close to the approximate thermal solution. In contrast, 
for the stable fixed point at late Langevin-time the plaquette averages are 
not time translation invariant along the contour — at least for the finite 
Langevin-times for which we followed the evolution. Here we employed a 
contour with r + = 0.5 and r_ = 3.5 with a t /a s = 0.03 and N t = 20 from an 
average over 10 runs. The real-time evolution near the thermal solution has 
been averaged over Langevin time 0.75 < $ < 1.25, for the non-unitary fixed 
point it has been averaged over 3.75 < $ < 4.25. 

As is shown is Sec. 14.21 all converging solutions of the stochastic 
dynamics method fulfil the same infinite set of (symmetrized) Dyson- 



31 




Figure 15: Numerical check of the Dyson-Schwinger equation for a spatial 
plaquette variable. Displayed are separately the LHS and the RHS of (|57J). 

Schwinger identities of the quantum field theory. This is remarkable in view 
of the different "physical" and "unphysical" solutions that are observed. In 
Fig-EHthis is visualized for the example of the Dyson-Schwinger equation for 
a spatial plaquette variable. Plotted are separately the LHS and the RHS of 
the Dyson-Schwinger equation (|57j) for the plaquette (see also Fig. [TJ. The 
plot displays the respective LHS and RHS as a function of Langevin-time. 
The flow with Langevin-time quickly leads to a rather accurate agreement of 
both sides such that the Dyson-Schwinger equation is fulfilled. However, after 
some Langevin-time they start deviating again, finally leading to another 
stationary value where the LHS and RHS agree to reasonable accuracy. Here 
the fluctuations of measured quantities of the system are much bigger. This 
is indicated by the given typical statistical error bars, with a comparably 
huge statistical error for the late Langevin-time evolution. We used median 
averaging for the right hand side of the equations. For the contour we employ 
r + = t_ = 2 with N t — 8 and a t /a s = 0.25. 
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7 Conclusions 



The motivation for this paper is the question of making real-time, 
nonequilibrium quantum field theory amenable to numerical simulations 
from first principles. This would not only boost our knowledge directly, 
it would also allow better testing of approximate analytical tools. This 
concerns, in particular, strongly interacting theories such as QCD, where 
reliable analytical approximations are difficult to find at phenomenologically 
relevant energies. Despite its importance research on real-time lattice gauge 
theory is still in its infancies. It has been a delicate problem and the 
various attempts based on reweighting or on Euclidean simulations have 
encountered major difficulties. Our aim here was to study the applicability of 
stochastic quantization to real-time, nonequilibrium problems. This method 
does, a priori, not involve any reweighting, nor redefinition of the Minkowski 
dynamics in terms of an associated Euclidean one. 

For setting up the procedure, both for scalar and for nonabelian gauge 
theory, we started from a five-dimensional classical Hamiltonian dynamics 
supplemented by a stochastic description and correctly accounting for the 
symmetries of the models. We then showed that the fixed points of the 
Langevin (fifth-time) dynamics fulfil the infinite set of Dyson- Schwinger 
equations associated with the respective quantum theory. We established 
in this way a set of identities to be fulfilled by the expectation values. These 
general results both allow checks of the simulations and suggest means to 
understand convergence and metastability problems. 

In the second part we undertook a numerical study of the real-time 
stochastic quantization approach applied to scalar and SU (2) gauge theory 
as a paradigmatic Yang Mills model. The main aim being to understand 
the problems of, and to develop means to control the method. Beyond 
varying the parameters of the implementation we also worked on various 
realizations of nonequilibrium and of non-zero temperature problems, which 
correspondingly define various integration contours. Our findings can 
be summarized as follows (both for scalar and Yang Mills fields, unless 
explicitely distinguished): 

- Instabilities of the Langevin dynamics are controllable: if the Langevin step 
size is chosen small enough run-away trajectories are rather seldom and the 
results do not suffer from discarding them. 

- Tilting the integration contour in the complex plane is a gauge invariant way 
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to improve both the convergence and the accuracy of the results. The physical 
effect of this tilt can be understood from exact results (e.g., Schrodinger 
equation for the anharmonic oscillator) and small tilts can therefore be used 
in the simulation. 

- Short real-time physics in thermal equilibrium can be reproduced reasonably 
well if the length of the real time contour is small on the scale of the inverse 
temperature f3, the details depending also on other parameters as explained in 
the main text. The Langevin flow is here dominated by a thermal fixed point 
which fulfils the unsymmetrized Dyson-Schwinger identities and describes a 
physical solution. This fixed point is stable for the scalar theory while it is 
metastable for the gauge theory, its life-time depending then on the contour 
and the other parameters of the problem. 

- For longer contours the boundary conditions in physical time do not seem 
to constrain enough the Langevin flow and the life-time of the thermal 
(physical) fixed point decreases (for gauge theory), or the fixed point 
becomes fully unstable (for scalar theory). A second, apparently stable fixed 
point develops for large Langevin times. This latter fixed point represents 
an unphysical, non-unitary regime, to be recognized by non-translational 
invariant expectation functions and violation of the unsymmetrized Dyson- 
Schwinger identities (while the symmetrized ones are still satisfied, indicating 
convergence). 

- For nonequilibrium problems similar observations hold. For the scalar 
model all these findings have been checked by comparing with results from 
a numerical solution of the corresponding Schrodinger equation. 

In view of these results we conclude that the Langevin approach can be 
used in direct numerical simulations of Minkowski theories, including Yang 
Mills theory, if a number of (rather severe) restrictions are observed. In 
those situations where physical solutions become unstable other, unphysical 
solutions will develop, but tests for differentiating them can typically be 
devised and applied. We did not try to optimize so far. The method allows 
for quite some flexibility as to which quantities are chosen to define the 
stochastic process, or introducing a stochastic reweighting. This is work in 
progress. 

Part of the numerical calculations have been performed on HELICS of 
the IWR Heidelberg. 
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